Entanglement Measures for Intermediate Separability of Quantum States 
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We present a family of entanglement measures Rm which act as indicators for separability of 
n-qubit quantum states into m subsystems for arbitrary 2 < m < n. The measure Rm vanishes 
if the state is separable into m subsystems, and for m = n it gives the Meyer- Wallach measure 
while for m = 2 it reduces, in effect, to the one introduced recently by Love et al. The measures 
R-m are evaluated explicitly for the GHZ state and the W state (and its modifications, the Wfc or 
Dicke states) to show that these globally entangled states exhibit rather distinct behaviors under 
the measures, indicating the utility of the measures Rm for characterizing globally entangled states 
as well. 

PACS numbers: 03.65.Ta, 03.65.Ud, 03.67.Mn 



I. INTRODUCTION 

Quantum entanglement typifies one of the most strik- 
ing aspects of quantum mechanics, posing profoimd ques- 
tions on our commonsensical comprehension of the phys- 
ical world. The conceptual significance of entanglement 
was first pointed out in the celebrated EPR paper 
where the nonlocal correlation of entangled states was 
regarded as a major obstacle for quantum mechanics to 
be a complete, realistic theory. The validity of nonlocal 
reality was later examined by Bell 0,01, who put the con- 
ceptual problem to one which is testable in laboratory. 
Since then, a variety of experiments have been conducted 
and by now we are almost convinced that non- 
locality does occur precisely as prescribed by quantum 
mechanics. Although the nonlocal correlation generated 
by quantum entanglement cannot be used for communi- 
cation 0, § I it suggests the existence of some nonlocal 
'influence' exerted between distant partners at a speed 
possibly exceeding that of light as reported by a recent 
experiment 0]. 

In view of its salient characteristics, quantum entan- 
glement is expected to play a vital role in our future tech- 
nology such as quantum computation and cryptography 
[13] ■ Successful application of entanglement will in gen- 
eral require the ability of manipulating and measuring 
entangled n-qubit states at a reasonable level of accu- 
racy. Among them, characterization of entanglement is 
perhaps the most basic requisite, and for this there have 
been a number of attempts including the use of canonical 
forms, entanglement witnesses and entanglement mea- 
sures [ill, [l3- These tools are certainly convenient for 
quantifying entanglement for a few small n cases, but 
they become almost intractable for large n due to the 
exponential increase in the number of distinct structures 
allowed for the entangled states [l^. It seems, there- 
fore, inevitable that in order to quantify entanglement of 
generic n-qubit systems, we need to resort to some means 
specifically designed for the objectives to be achieved. 

Among the many entanglement measures proposed so 
far [ll|,[li,[i3,[il,[li, the Meyer- Wallach (MW) measure 
[r^ is notable in that it examines the full separability, 
i.e., if the n-qubit state under inspection is a product 



state of all the n constituent subsystems [3 [3- Re- 
cently, Love et al. [2^ proposed a measure which is 'oppo- 
site' to the MW measure in the sense that it examines the 
global entanglement, i.e., if the state admits no two sub- 
systems into which it can be decomposed as a product. 
In the present paper, we present a family of entanglement 
measures Rm, rn = 2,3, ... ,n which can examine the in- 
termediate separability, i.e., if the state is a product state 
of arbitrary m subsystems. In particular, for m = n our 
measure coincides with the MW measure, whereas for 
771 = 2 it reduces, in effect, to the measure of [1^. We 
show that, besides as indicators of intermediate separa- 
bility, our measures can also be used in characterizing 
globally entangled states in general. This is illustrated 
by the two standard globally entangled states, the GHZ 
state 21 1 and the W state [ij] in n-qubit systems, which 
exhibit rather contrasting behaviors under our measures 
Rm for various m. Analogous distinct behaviors can also 
be observed for the set of globally entangled Wfe states 
(known as Dicke states), which are introduced as modi- 
fied W states for 1 < k < [n/2j, with W„/2 furnishing 
the maximally entangled state for the MW measure i?„. 

The present paper is organized as follows. In the next 
section, we introduce the family of entanglement mea- 
sures with the required intermediate separability, and 
show that both the MW measure and the measure of 
[20| appear at the two ends of the set. We then analyze, 
in section 3, the globally entangled GHZ and W states in 
terms of the measures introduced. In section 4, the anal- 
ysis is extended to the Wfe states. Section 5 is devoted 
to our conclusion and discussions. 



II. ENTANGLEMENT MEASURES AS 
INDICATORS OF INTERMEDIATE 
SEPARABILITY 

The system we consider is an n-qubit system whose 
quantum states are described by vectors in the Hilbert 
space . In order to discuss its arbitrary subsystems, 
we label the n constituent 1-qubit systems by integers so 
that any subsystem consisting of some of the constituent 
systems is specified by a subset of T = {1, 2, . . . , n}. Let 
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V = {si\^i be a partition of T, 



2.6. 
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T and s, n 



for i 7^ J. 



(1) 



Each subset Si determines a corresponding subsystem of 
the total system , and hence we may use Si to refer 
to the subsystem specified by the subset. We denote by 
Si the subset complementary to Si in T with Si Usi~T. 

Now, given a pure state let be the reduced 
density matrix in the subsystem Si obtained by taking 
the trace of the density matrix p = over the com- 

plementary space Si. Letting also \si\ be the number of 
elements (constituents) in the subset Si, we recall that 
the quantity 



7jJtlj) = N{\s,\){l-typl), N{\s,\) 



- 1' 



(2) 



introduced in [20] vanishes risiO^) = iff the state {tp) is 
separable with respect to Si and s^. Here, the normal- 
ization factor iV(|si|) in ^ is chosen so that we have 
Vsi (V') 1 when the reduced state is maximally mixed 
Psi — j^-^- This quantity rjs- is in fact a generaliza- 
tion of the (squared) concurrence [l^l and can also be 
regarded as the quantum linear entropy [23j . From the 
quantities rig. in ([2]) obtained for all the subsets Si in V, 
we evaluate the 'average' value for the partition V by the 
arithmetic mean. 



(3) 



Clearly, we have ^p(^) = iff the state jV') is separable 
according exactly to the partition V of the total set T. 

Out of all possible partitions V of T, we may choose 
those V consisting of m subsets for some m in the range 
2 < m < n, and evaluate the geometric mean of the quan- 
tities ^-p (?/;). Namely, if (i(7^) is the number of subsets of 
the partition V ^ we consider 



1/S(ri,m) 



-R™(V') 



n ^^(^) 



where 



S{n,m) = ^ 



-1) 



m — k h.n — 1 



^ (fc- l)!(m-fc)! 



(4) 



(5) 



is the Stirling number in the second kind |24| , which rep- 
resents the number of all possible partitions of the inte- 
ger n into m subsets, or the number of partitions V with 
(1(1^) = m. The quantities Rm{ip) possess the important 
property: 



i?™(^) =0 ^ 



1^) is separable (at least) 
in TO subsystems in . 



(6) 



Besides, since i?m('0) ^-re formed from r]si{tp) which are 
all entanglement monotones [l^l , each of them, Rm {tp) , 
TO = 2, . . . , n, qualifies as an entanglement measure. In 
particular, for m = n where the subsets Si, i = 1, . . . , n, 
correspond to all the constituent subsystems, we have 



n 

RnW = -Y.2{l-trpl}, 



(7) 



which is precisely the MW measure [13, [13 (see also [2^, 
[27|V On the other hand, at the other end to = 2 we have 
the partitions S = {si = s, S2 = s}. Choosing the subset 
s so that \s\ < \s\, and noting S{n, 2) = 2"~^ — 1, we find 



1/(2—1-1) 



R2W=\ n c(s)77.(^) 

Ki<\s\<\s\ 



(8) 



where the prime on the product symbol indicates that 
either one of the subsets s and s is included when \s\ = |s|, 
and the coefficients c(s) are given by 



1 2l''l - 2!"l 

c s =1 ^. 

w 2 2" - 2l^l 



(9) 



The measure R^bP) is equivalent to the measure proposed 
by Love et al. [20|, apart from the factor c(s) which varies 
between the maximum c(s) = 1 for \s\ ~ \s\ and the 
minimum c(s) = 1/2 -|- 1/(1 - 21"") > 1/2 for \s\ ^ 1. 

For illustration, we consider, for example, the n-qubit 
state |Gfc)|0)®("~'^) with |Gfc) being the fc-qubit version 
of the GHZ state 

|GHZ) = i= (111... 1) + 100...0)). (10) 

The state |Gfc)|0)®("~''') is clearly separable into n— fc-f-l 
subsystems, and Figure[T]shows that the measures Rm{ip) 
for fc = 3, . . . 8 and n — 8 indeed vanish at to = 8 — fc + 1 , 
where we further observe that i?j„ behave rather dis- 
tinctively depending on the size k of the entangled sub- 
system. The measures are also evaluated for the prod- 
uct states |Gfc)|G/) consisting of two GHZ states for 
{k, I) = (4, 4), (5, 3), (6, 2), and the result shows that Rm 
can distinguish these states which are all indistinguish- 
able under both of the MW measure (since Rg = I) and 
the measure by Love et al. (since i?2 = 0) . These observa- 
tions suggest that our measures Rm{ili), as a whole, may 
also be useful to characterize multipartite entanglement 
of the state IV") in addition to examining the intermediate 
separability. This possibility is explored further in terms 
of the GHZ state and the W state later, where we also 
present an explicit procedure to evaluate the measures 
for these particular states. 

At this point, we mention that the quantity rjs. in p|). 
which is an entanglement measure for the separability 
of the subset s^, can be considered as the purity mea- 
sure based on the subalgebra associated with Sj in the 
generalized framework of entanglement introduced ear- 
lier [2^, [2^, [231 ■ It is also worth mentioning that the 
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FIG. 1: The values of the measures Rm evaluated for two 
sets of 8-qubit states: |Gfe) |0)®*®"*' for fc = 3, . . . , 8 (above) 
and |Gfc)|GO with {k,l) = (4, 4), (5, 3), (6, 2) (below), which 
are represented by the tags 'Gfc' and 'GkGV, respectively. 
The vanishing values confirm the number of subsystems into 
which the states |Gfc)|0)®*""''' are factorized. AU of the states 
\Gk)\Gl) which are separable only into two subsystems share 
the same values ii2 = and iig = 1 but allow different values 
for other Rm- 



same quantity r]s- , or more generally the mean value ^-p 
in ([3|) evaluated for the given partition P, is related to the 
quantum Fisher information for the parameter estima- 
tion of the low-noise locally depolarizing channels whose 
actions for quantum states are specified by the partition 
V. Since the inverse of the quantum Fisher information 
gives the lower bound of the variance of estimators, we 
can provide the operational meaning for ^-p as a mea- 
sure of precision in the estimation of the strength of the 
low-noise locally depolarizing channels associated with V 
[2^. This relation between and the quantum Fisher 
information implies that Rm may be interpreted as the 
quantum Fisher information for an assembly of the low- 
noise depolarizing channels under the condition that only 
the number of the local channels is known. 

The computational complexity of Rm may be esti- 
mated based on the simple rule that all arithmetic oper- 
ations (addition, multiplication, division and taking the 
fc-th root) are equally counted. We then find that, since 
the number of summations needed for ps- is 2"~l*'l, the 
computational complexity of Rm grows exponentially in 
general. Note that this applies to any measures (includ- 
ing the MW measure) which require the partial trace 
operations. Further, if rjs-ijp) for all Si are given, the 
number of additional steps necessary to obtain Rmi^p) is 
approximately 'mS{n,m). Since S{n,m) grows exponen- 
tially for large n with fixed m (see the Appendix) , so does 
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FIG. 2: The schematic diagram of scalable sequences of Rm 
for computational complexity. The dot on the lattice point 
(m, n) represents Rm of n-qubit states. Dotted arrows indi- 
cate the polynomial growth required to obtain Rm from the 
given set of rja^ , while solid arrows indicate the exponential 
growth. The MW measures R„ lie on the lower diagonal edge 
of the triangle area of the lattice while the measures R2 lie on 
the vertical left edge. 



the number of steps for Rm{ip) from T]s^{ip). On the other 
hand, if we are interested in the measures i?„_; {4^) with 
fixed I for large n {e.g., the MW measure Rni'ip) arises 
for / = 0), then we see that the required number of steps 
is 0(71,^'+^), that is, it grows only polynomially. These 
observations show that the computational complexity of 
Rmi'4') for 'T'-qubit states depends on how to construct 
the scalable sequences of the measures in question (see 
Figure [2]). The polynomial growth will also arise in gen- 
eral when we restrict ourselves to symmetric states [29|. 

The measures Rm{ip) can readily be extended to those 
which accommodate mixed states as well [s^. This is 
done by adopting the standard procedure of considering 
the convex hull at the stage of ^-p : 



iv{p) = ^ min 'Vpa^p(V'a), 



(11) 



where the minimum is chosen from all possible decompo- 
sitions of the density matrix p = ("0111 into the 
probability distribution {pa} and the pure states |0a). 
From the extended S,v[p) in pT|) . we define Rm[p) as 
(III). The resultant measures Rm{p) possess the desired 
property of intermediate separability as an extension of 
^ to mixed states. Namely, Rm{p) = iff the state p 
is separable into m subsystems as a mixed state, i.e., it 
admits the form, 



(12) 



i=i 



where p". are density matrices in the subsystems Si . Note 
that so defined Rm{p) become entanglement measures, 
since Rm{p) are monotone for LOCC and invariant under 
local unitary operations. 
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III. GHZ STATE VS W STATE 

We now evaluate the amount of entanglement pos- 
sessed by the two familiar globally entangled states, the 
GHZ and the W states, using the measures Rm{ip) in- 
troduced above. These are particular states which are 
invariant under all permutations of constituent subsys- 
tems, and this exchange symmetry facilitates our com- 
putation considerably. To proceed, we first note that for 
those symmetric states the quantity 77s. in ([2]) depends 
only on the number of the elements \si\ of the subset 
Si, not on the choice of the elements in s^. To find the 
value of the measure Rm{ip), we need to consider all pos- 
sible partitions V with = m to get the quantity 
S,v{ip) in Q, but again the exchange symmetry implies 
that ^-p('0) depends only on the way the partition V is 
formed in terms of the set of numbers |sj;| of the elements 
in the subsets Si comprising V . To be more explicit, let 
us choose the numbering of the subsets s,; in the order 
< |s2| < • • • < \sm\ and introduce the notation, 



^\}■ 



(13) 



Note that \V \ furnishes an ordered partition of the integer 
n into d{V) = m nonvanishing integers by n = |si| -I- 
|s2| + • • • + |sin|- For n and m with 2 < m < n, let 
0{n,m) be the set of all distinct ordered partitions of 
the integer n into m nonvanishing integers. Given some 
1^1 € Q{n,m), we denote by h{\V\) the total number 
of partitions V sharing the same ordered partition |7^|. 
The measure in O can then be calculated by the 

product of ^^(-0) for all different \V\ in Q{n,m), i.e.. 



1 / S{n.rn) 



n i^wim 



H\v\] 



(14) 



where we have written C|7'|('0) for ^vii') to stress that it 
is dependent only on \V\. 

Now we consider the n-qubit GHZ state (fTO|) . The 
GHZ state is quite special since it has tr(psj^ = 1/2 for 
all subsystems s^, and from this we obtain 



77,, (GHZ) = 



(15) 



with A^(|si|) given in ([2]). To illustrate our procedure 
for evaluating the measures, we choose, for instance, the 
case n = 4, TO = 2 for which the set Q{4, 2) consists of the 
two elements, \V\ = {1,3} and {2,2}. The numbers of 
partitions with the same \V\ are, respectively, h{{l, 3}) = 
4!/(l!3!) = 4 and h{{2,2}) = 4!/(2!2!2!) = 3, yielding 
5(4,2) = E|P|ee(4,2) ^(1^1) = 7. We then find 



i?2(GHZ) = - 




0.732. 



(16) 



This procedure can be applied for any n and to, and the 
results up to n = 50 arc shown in Figure [31 
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FIG. 3: The entanglement measures Rm evaluated for the 
GHZ state |GHZ> (above) and the W state |W) (below) as 
functions of m for various n with 3 < n < 50. Each curve 
represents Rm for m in the range 2 < m < n with a fixed 
value of n which can be read off from the right end value of 
m of the curve. 



Next, we consider the W state. 



|W) = ^ (|10 • • • 0) + |01 • • • 0) + • • • + |00 • • • 1)) , (17) 



which has 



,7,.(W) = 7V(|s,|) 



2|s.|(n- |s,|) 



(18) 



For comparison, we again choose the case n = 4. to = 2 
to find 



i^2(W)= - 



33\^ /2-31^ 



0.621, 



(19) 



which is less than the value of the GHZ state. As 
in the GHZ case, the results up to n = 50 are shown in 
Figure [3l 

It is clear from Figure [3] that the GHZ and the W 
states exhibit rather contrasting behaviors for the en- 
tanglement measures R„i. Namely, for the GHZ state, 
Rm is a monotonically increasing function of m confined 
within 1/2 < Rm < 1 and approaches the value i?„ = 1 
at the right end m = n. In contrast, for the W state, 
Rm is basically a decreasing function of to confined in 
< Rm < 1/2, except for the small n < 9 for which 
Rm can exceed the value 1/2. These can also be seen 
directly from the formulae (|15|) and ^TE\\ . In a sense. 
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this result agrees with our intuitive picture of the GHZ 
state being more globally entangled than the W state for 
all n, which is also observed by using the entropy of en- 
tanglement [2^. Another point to be noted here is that 
the clear difference in the values of Rm between the two 
sets of states suggests that the GHZ state is more fragile 
than the W state, because the measures are indirectly 
related to the fragility of the state which is the source of 
the operational meaning of preciseness in the estimation 
discussed in [2^. This again is consistent with the con- 
ventional view that the entanglement of the W state is 
more robust than that of the GHZ state [10, [3l[ . This 
propensity of robustness may also be recognized by con- 
sidering relevant combinations of partitions V specific to 
that purpose as done in poj . 

Note that the lower bound 1/2 of the measures Rm 
for the GHZ state indicates that the GHZ state cannot 
be approximated well by a state which is separable in m 
subsystems for any number of m. On the other hand, for 
the W state we observe that the values of Rm with m 
closer to n approach zero for larger n, and in particular, 
the value i?„ (i.e., the MW measure) has the vanishing 
limit, 



Rm(Wio) 



lim Rn{W) = lim 



1 



0. 



(20) 



This, however, does not mean that the W state |W) be- 
comes fully separable in the limit n ^ 00. One can see 
this by considering a geometric measure of entanglement 
Eq{iI)) which is defined by 



Sg(^) := l-max|(x|V)l^ 

X 



(21) 



where the maximum is taken over the set of fully separa- 
ble states Ix) [13, m, [13, [13 ■ Indeed, parameterizing an 
arbitrary fully separable n-qubit state as 



Ix) =(8)(cos0,|O), + e'^' sin0,|l),). 



(22) 



and varying the angle parameters in |x), one finds that 
the value of Eq ( W) is obtained when sin^ 9i = 1/n for 
all i and 4>i = for all Hence, in the large n limit 
we find 1351 



lim Eg{^) = 1 



lim ( 1 

n — ^00 



1 



(23) 



which shows an intriguing fact that despite the vanishing 
limit of the MW measure i?„(W), the W state does not 
approach a definite fully separable state in the limit n — > 
00. This indicates that the connection (|4]) between the 
vanishing measure and the separability, which is perfectly 
valid for finite n, does not hold for n 00. 



IV. MODIFIED W STATES 

We may further examine the property of our measures 
by considering a set of states which arc totally symmetric 
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FIG. 4: The entanglement measures Rm evaluated for the 
Wio state (above) and the W20 state (below) as functions of 
m for various n with fc -I- 1 < n < 50. Each curve represents 
Rm for m in the range 2 < m < n with a fixed value of n 
which can be read off from the right end value of m of the 
curve. 



with more than one |1) states in the constituent subsys- 
tems. To be explicit, we introduce the 'Wfe states' for 
l<k< [n/2\ by 



|Wfc 




perm. 



(24) 



where 'perm.' means that all possible distinct terms pos- 
sessing fc 'I's and {n — k) 'O's obtained by permutations 
of the first term are included. The Wfc states, which are 
known as Dicke states, reduce to the standard W state 
|Wi) ~ |W) for fc = 1, while for fc > 1 they become 
slightly more involved but are still manageable thanks to 
the symmetry. 

To evaluate the measures, we first implement an ap- 
propriate unitary transformations to |Wfc) so that, for 
a given subsystem Si, the state in Si is represented by 
the left I Si I qubits in the n-qubit state |*) = |*)sj*)si. 
To proceed, it is also convenient to specify each of the 
terms in |Wfc) by the number of 'I's, which is fc for 
|Wfc), and an integer ct for 1 < tr < ('^) labeling the 
distinct terms appearing in the permutations. Clearly, 
the same notation can be employed for both of the sub- 
systems Si and Si as well, and we may write an arbitrary 
term in (|24p as a product of states in the two subsys- 
tems as Ifc,^) = |r, r)sjfc — r, t')s; , where |7',r)s. is a 
state of the subsystem with r 'I's and the label t 
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FIG. 5: (Color online) The entanglement measures Rm eval- 
uated for the Wfe states as functions of m and k under the 
fixed value of n = 40. 



runs for 1 < r < ('^''), and similarly \k — r, t')^; is a 
state of the subsystem Si with the label r' running over 
1 < t' < ("^ri^''). This allows us to rewrite the state 
(l24ll in the form, 



|Wfe> 



(25) 



from which the reduced density matrix is found as 

"EfTJ!'')El^' -'I' (26) 



where the summation of r is for max(|si| — (n — fc), 0) < 
r < min(|si|, k). It is now straightforward to evaluate rj^. 
to find 



E 



.(27) 



Based on the result (j27p . one can obtain the values of 
the measures Rm for the states, and the outcomes are 
shown in Figure [4] for the two cases A: = 10 and k ~ 20. 
It is seen that both of the Wio and W20 states exhibit 
distinctive behaviors which are also different from those 
of the GHZ states and the W states discussed before, and 
in particular we notice that the measure Rm achieves the 
upper limit Rm = 1 at n = m = 2k. This can also be 
confirmed from (j27p since for symmetric states, n = m 
implies Rm{'4') = VsAi^) with \si\ = 1, which is 1 for 
n = 2k. We therefore see that the n-qubit W„/2 state 
furnishes the maximally entangled state for the MW mea- 
sure. This may be understood from the large symmetry 
possessed by the W„/2 state which has an equal number 
of |0) and |1) states in the constituent subsystems. 

It is also interesting to look at the behaviors of the 
Rm measures for the Wfc states as functions of k with 
some fixed n. This can be done in Figure [SJ where we 
plot the values of Rm for Wk in the range 2 < to < 40 
and 1 < A: < 20 for n = 40. We observe there that 
Rm arc monotonically increasing functions of k for all 



TO, indicating that the Wk states are 'more entangled' 
for larger k < n/2 under all measures Rm- Moreover, 
we see that the change in the values of the measure is 
in general more prominent for Rm with higher to, which 
suggests that variation of k alters the Wfc states in their 
entanglement property of higher separability. 



V. CONCLUSIONS AND DISCUSSIONS 



In this paper we have introduced a family of entangle- 



ment measures R„ 



for n-qubit states in 



which both the MW measure [l^ and the measure pro- 
posed by Love et al. [l^l arise as the two extreme cases 
m = n and m = 2. Our measures Rm are scalable and 
can be used as indicators for separability of the n-qubit 
states into to subsystems. We have seen by comparing 
the behaviors of the GHZ state and the W state that the 
measures Rm arc also useful to characterize the entan- 
gled nature even for globally entangled states. The Wfe 
states, which give the maximally entangled state for the 
MW measure at k = n/2, furnish another set of globally 
entangled states which behave distinctively under Rm de- 
pending on k. 

The outcomes of the analysis on the GHZ state and 
the W state indicate that, in terms of our measures Rm, 
the GHZ state is more entangled than the W state. This 
is to be contrasted to the observation in [H, [111 that, 
in terms of the geometric measures, the converse holds. 
This may be derived from the difference in the basic in- 
gredient of the measures, that is, for bipartite systems our 
measures Rm are a generalization of the linear entropy, 
whereas the geometric measures are related to Chebyshev 
entropy [2^. Apart from the origin of the difference, this 
poses the question on the relation between the two fam- 
ilies of entanglement measures, and calls for clarification 
of the physical and operational meanings of these mea- 
sures other than the aforementioned preciseness of esti- 
mation valid for the constituent element ^p. Along with 
the formal classification of our measures with resp ect to 
the generalized framework of measures [H, [2^ |27| |. these 
issues should be investigated further to gain a fuller pic- 
ture of multipartite quantum entanglement in general. 
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APPENDIX 

This appendix provides an elementary account of the 
asymptotic behavior of the Stirling number in the second 
kind S'(n, to) for large n, which is needed to estimate the 
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computational complexity of Rm{ip) in the text. Recall 
first that S{n, m) fulfills the recurrence relation, 

S{n, m) — S{n — l,m — 1) ~ mS{n — f , to) (A.f ) 

with the initial conditions S{n,n) = S{n, 1) = I. From 
this we immediately obtain 

S{n, to) > mS{n -l,m)> ■■■> to""", (A.2) 

which shows that, for fixed to, the number of steps to 
compute Rm{ip) grows exponentially as 0{m"-). 

Next, we consider the case where n increases with fixed 
I as discussed in Sec. II. For example, if Z = 2, we sum 
up the recurrence relations (jA.ip with n~k,in — k~2 
from fc = 4 to fc = 71 and use S{n, n — I) = (^j to obtain 

S{n, n - 2) = ^n{n ~ \)(n ~ 2)(3n - 5) = 0{n^). (A.3) 

More generally, if we assume S{n,n ~ (l ~ 1)) = 
©(n^^'"^)), then by the same iterative procedure for 



()A.ip . we obtain S{n,n — I) = 0{n'^^). By induction, 
this gives the polynomial growth of S{n,n — I) for large 
n. 

With the help of Stirling's approximation, a more 
explicit asymptotic formulas of S{n, to) can be found. 
Namely, for m, S{n, to) is written as [sH] 

which shows the exponential growth of S{n, to) for large 
n and is consistent with our result S{n,m) = ©(to"). In 
[st! ]. we also find 



S{n,m)^ ' (A.5) 
[n — my. 

for large n and small n — m, which confirms our result 
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